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Acceleration from Modified Gravity: Lessons from Worked Examples 
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University of Chicago, Chicago, IL 60637 

I examine how two specific examples of modified gravity explanations of cosmic acceleration help us understand some 
general problems confronting cosmological tests of gravity: how do we distinguish modified gravity from dark energy if they 
can be made formally equivalent? how do we parameterize deviations according to physical principles with sufficient generality, 
yet focus cosmological tests into areas that complement our existing knowledge of gravity? how do we treat the dynamics 
of modifications which necessarily involve non-linearities that preclude superposition of forces? The modified action f(R) 
and DGP braneworld models provide insight on these question as fully-worked examples whose expansion history, linear 
perturbation theory, and most recently, non-linear A^-body and force-modification field dynamics of cosmological simulations 
are available for study. 



1. INTRODUCTION 

The next generation of surveys aimed at characteriz- 
ing the acceleration of the expansion will also provide 
precision tests of gravity in a cosmological context. In 
particular, they will address the possibility that the ac- 
celeration itself is driven not by an unknown dark energy 
component but rather by a modification to gravity on 
the largest scales (e.g. |H2l3l4j ). 

In the absence of a compelling theory for modi- 
fied gravity, how do we best quantify and focus these 
tests? Is there an analogue of the parameterized post- 
Newtonian description of local gravity for cosmology? 
Can it act as a meeting point between observations and 
theory? 

At the very outset there are some basic questions 
that confront this endeavor. In particular, three issues, 
loosely based on the challenges compiled in 0], appear 
as obstacles: 

1. Dark energy equivalence. Formally, any met- 
ric theory of modified gravity can be recast into 
an equivalent dark energy theory under ordinary 
gravity. Is a cosmological test of gravity itself even 
possible? 

2. Anything goes. Without the framework of general 
relativity we are left with arbitrary functions of 
space and time to relate the matter and the metric. 
What guidance do we have to parameterize them? 
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3. Superposition. Order unity cosmological devia- 
tions and tiny local deviations imply that the mod- 
ified Poisson equation is non-linear in some way. 
How do we treat non-linear dynamics without a 
superposition principle for gravitational forces? 

Given such basic issues, it is useful to have concrete 
examples, however contrived, ad hoc, or problematic on 
a quantum level. They serve as case studies for how a 
modified gravity theory might work. 

Two such examples are the modified action model 
where the Einstein-Hilbert action is supplemented by 
a non- linear function of the Ricci scalar f(R) [5161 
[7] and the 5-dimensional braneworld model of Dvali, 
Gabadadze, & Porrati (DGP) [8]. Both cases are now 
fully worked in that their expansion history, linear per- 
turbation theory, and most recently, non-linear TV-body 
dynamics of cosmological simulations have been studied. 

In this piece, I examine how these worked examples 
provide insight on the three conundrums confronting 
cosmological tests of gravity for cosmic acceleration. 

2. DARK ENERGY EQUIVALENCE 

Given the formal equivalence which maps a modified 
gravity theory onto a dark energy theory under ordinary 
gravity, can cosmological observables ever test gravity 
itself? 

Let's examine the formal equivalence and then place 
it in the context of the two worked examples. A metric 
theory of gravity where energy- momentum is covariantly 
(locally) conserved (V M T Ml/ = 0) can always be rewritten 
as a dark energy theory. To see this, consider that a 
generally covariant field equation of the form (e.g. [9ll0j ) 
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can be re-expressed as ordinary gravity G^ — SttGT^ 
with an effective dark energy stress energy 
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Moreover this effective stress-energy is covariantly con- 
served V M T^ E = by virtue of the Bianchi identity. 
Also, conservation of the matter means that the only 
influence of the modification to gravity or effective dark 
energy comes through the metric: once the gravitational 
potentials are determined, matter falls in them as usual. 

While this formal equivalence holds, it does not mean 
that the equivalent dark energy description obeys rea- 
sonable microphysics. Microphysics determine the equa- 
tions of state that relate the components of the stress 
energy tensor: the pressure, energy density, momentum 
density and anisotropic stress. The effective dark en- 
ergy description of a modified gravity model will look 
contrived as its microphysics must mimic an explicit de- 
pendence on the metric unlike in physical dark energy 
models like scalar fields [ID] . 

Moreover, even though the matter and effective dark 
energy are separately conserved, the relationship be- 
tween the effective microphysics and the metric imply 
an implicit coupling between the dark energy and mat- 
ter - i.e. the effective dark energy mediates a fifth forces 
that universally couples to the matter. The f{R) and 
DGP models arc both exhibit this phenomenology. They 
contain an extra scalar degree of freedom that can ei- 
ther be represented as dark energy anisotropic stress, 
proportional to the metric, or the propagating mediator 
of a fifth force. As such, they are phenomenologically 
quite distinct from typical dark energy models. 

In summary, cosmological tests of gravity, just like so- 
lar system measurements, require a notion of what com- 
poses the stress-energy tensor (s) in the system. For cos- 
mological tests, we assume candidates for the dark mat- 
ter and dark energy with microphysics that, along with 
conservation of energy-momentum, closes their equa- 
tions of motion [IT] . What most cosmologists implicitly 
mean by a "test of gravity" is: a test of gravity in the 
context of cold dark matter and dark energy whose den- 
sity is relatively spatially smooth on scales well below 
the current horizon^ In this context, general relativity 
makes specific and falsifiable predictions for cosmolog- 
ical observables. Falsification however does not neces- 
sarily imply that general relativity is incorrect. 

3. ANYTHING GOES 

Without the framework of general relativity how do 
we parameterize and test gravity cosmologically? 



Consistency Approaches: If anything goes, then 
shouldn't one start with parameters that define the 
simplest self-consistency tests of general relativity (plus 
dark energy that is smooth relative to cold dark mat- 
ter)? 

Candidates for self-consistency tests include an ap- 
proximate expression for the growth rate of linear den- 
sity perturbations A m given the matter contribution to 
the expansion rate Q, m (a) (e.g. [12113] ) 
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2 Modified gravity explanations of the dark matter fall into a dif- 
ferent category and are not the subject of this exploration. 



where 7 « 0.55. Other approaches include separate dark 
energy equation of states that govern growth and dis- 
tance measures (e.g. [14 15 ) or model independent re- 
constructions of growth from distance measures that can 
be compared with growth measurements (e.g. |16|17j ). 
The latter have the advantage that consistency checks 
can be made exact in linear theory whereas the former 
two parametric schemes are easier to apply. 

One problem with these approaches is that deviations 
in these parameters away from consistency do not rep- 
resent physically possible modified gravity models. As 
minimal approaches, they also do not target areas where 
measurements might best complement what we already 
know about gravity. For example, consistency statistics 
can be dominated by information from the non-linear 
regime where, as we shall see, ordinary gravity is al- 
most guaranteed. Conversely, an incomplete account of 
non-linear baryonic physics can give false consistency 
violations [18] , 

Likewise, characteristic signatures of modified gravity 
aren't very well exposed by these consistency statistics. 
Neither the f(R) nor DGP models have a scale indepen- 
dent linear growth rate, though deviations for the latter 
only appear near the horizon scale [19120121] , Moreover 
modified gravity theories such as these change the re- 
lationship between the gravitational potential, spatial 
curvature and the density perturbations. 

It is tempting to simply supplement consistency pa- 
rameters like 7, generalized perhaps to be a function of 
scale, with a variation in the Newton constant G that 
also depends on scale (e.g. [22]). However these types of 
prescriptions do not establish a framework where phys- 
ical principles like conservation of energy-momentum 
strictly apply nor do they provide a viable construction 
beyond linear perturbations. 

Despite these flaws, simple consistency tests are valu- 
able in that they are easy to apply and forecast. Their 
violation will point to new physics or astrophysics that 
can then be more fully developed. They are not however 
the only reasonable metric to apply for judging tests of 
gravity. 

Linear Parameterization: Fortunately, the problem of 
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formal equivalence to dark energy in $2] and knowledge 
that gravity must only have small deviations from gen- 
eral relativity locally (see Q is a virtue for developing 
a parameterized description of modified gravity. 

The framework of general relativity and aspects of 
the parameterized post Newtonian description for devi- 
ations carry over to modified gravity. Coupled with spe- 
cific worked examples of the types of deviations expected 
in modified gravity models, this framework brings form 
to a cosmological parameterization of gravity. 

As discussed in Sj2] if we can parameterize the degrees 
of freedom associated with the microphysics of the effec- 
tive dark energy then we have a description that satisfies 
physical principles: gravity as a metric theory and and 
conservation of energy-momentum. 

Modified gravity theories, including the two worked 
examples, typically have an extra scalar degree of free- 
dom associated with the physics that drives the accel- 
eration. Just as in the parameterized post-Newtonian 
(PPN) description, this scalar characterizes the rela- 
tionship between the Newtonian potential '5 and the 
spatial curvature $ but now in the context of small de- 
viations from a statistically homogeneous and isotropic 
Friedmann- Robertson- Walker background 



ds 2 = -(l + 2*)cft 2 + a 2 {l + 2$)dx 2 . (5) 

Note that g = for a universe composed of cold dark 
matter and smooth dark energy. These deviations can 
be measured by comparing gravitational lensing and 
gravitational redshifts, which are sensitive to <I> — "f, 
with the dynamics of non-relativistic matter, which is 
sensitive to ^. 

In terms of the effective dark energy, the metric devi- 
ation g is related to anisotropic stress (e.g. [33 24J). On 
superhorizon scales, this parameter g(a, k — > 0) alone is 
sufficient to describe the evolution of the two potentials 
\& given a known, separately parameterized, expan- 
sion history [25] . Completeness follows from the conser- 
vation of the comoving curvature for adiabatic pertur- 
bations 

C = const, (k < aH) (6) 

independently of the microphysics of dark energy. This 
independence can also be directly shown without the 
aid of the dark energy equivalence [26j . Fundamentally, 
it only requires that conservation of energy-momentum 
still apply. 

This fact is also quite useful for describing the two 
worked modified gravity examples. In the DGP sce- 
nario, the complicated superhorizon physics of met- 
ric fluctuations propagating into the extra dimension 



through the master equation for the bulk is simply en- 
capsulated by this one function on the brane [27] . The 
inclusion of this horizon scale effect leads to an excess 
in the large angle CMB anisotropy [35J which brings the 
self-accelerating solution into 5a tension with the joint 
CMB and supernovae data [29 . Likewise, the compli- 
cated 4th order dynamics of f(R) theories is also neatly 
encapsulated in g on large scales [T9"] . 

On smaller scales, conservation of energy-momentum 
and g is not sufficient for a complete description. Again 
dark energy equivalence is useful for understanding this 
fact. There it is necessary to close the equations of mo- 
tion by supplementing conservation with an equation of 
state that relates the pressure fluctuation to the den- 
sity, momentum and anisotropic stress fluctuations. For 
a smooth dark energy component, this is achieved by 
setting a Jeans scale or a sound speed for pressure and 
anisotropic stress fluctuations for the dark energy [llj . 

One possibility that has been explored |24l30j is set- 
ting the anisotropic stress and pressure of the effective 
dark energy to exactly cancel so as to leave no source 
to the momentum density in a particular gauge choice, 
synchronous gauge. This amounts to assuming that £ 
is a constant on all scales if the remaining matter is 
cold dark matter and corresponds to a case where the 
effective dark energy has non-negligible density contri- 
butions on small scalesH However it is too restrictive 
and does not lead to a small scale description that has 
the features of the worked examples. 

In both the DGP and viable f(R) models, the New- 
tonian regime is well approximated by an unmodified 
Poisson equation for the lensing potential 

V 2 f^-^) = -^Ga 2 p m A m , (7) 

where spatial derivatives are comoving throughout. 
Along with g and energy-momentum conservation for 
the matter, the Poisson equation for the lensing po- 
tential completes the system. More generally G can 
be made time dependent without breaking the struc- 
ture of these equations [10]. On the other hand, pa- 
rameterizations that just involve g and the Poisson 
equation 31 2 2132] fail to automatically enforce energy- 
momentum conservation or the Bianchi identity on hori- 
zon scales and above. For the f(R) model g — » —1/3 
below the Compton wavelength scale in the background 
and for the DGP model g — > — l/3/3(a) where /3(a) 
1 — 2Hr c (l + -din H/dlna) with r c as the crossover 
scale. 

3 The same issue applies to synchronous gauge techniques that 
simply shut off dark energy perturbations on small scales as an 
approximation to smooth dark energy: in a synchronous gauge 
choice where the dark matter defines the frame, the momentum 
density of the dark energy is never negligible. 
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A parameterized post-Friedmann (PPF) description 
that is flexible enough to incorporate both the large and 
small scale behavior of linear fluctuations in the two 
worked examples for modified gravity was introduced in 
|fOj . The gist of this parameterization is to include an 
interpolation between the large scale Friedmann behav- 
ior of equation ([6]) and the Newtonian behavior of equa- 
tion (|7|) while strictly conserving energy-momentum on 
all scales. This amounts to highly non-trivial closure re- 
lations or equations of state for the equivalent effective 
dark energy. Such complexity addresses the question of 
distinguishability raised in $2] 

These relations are readily generalized to the early 
universe where radiation fields from the photons 
and neutrinos are important and provide their own 
anisotropic stress [33]. A public code to calculate lin- 
ear fluctuations including the CMB under this general- 
ized PPF description was released in 34J This method 
has also been adapted to handle smooth dark energy, 
braneworld gravity with brane tension [35] . and explo- 
rations of degravitation ideas [36] . 

Non-linear Parameterization: There is one final, per- 
haps most critical, piece to a PPF description. Given 
that there exists very strong constraints on analogous 
PPN parameters locally, e.g. in solar system tests, there 
must be some means by which finite modifications at 
cosmological scales are reduced to acceptable levels lo- 
cally 

Here is where having the worked examples prove the 
most valuable. In the DGP and f(R) scenarios, ordi- 
nary gravity is restored by having non-linear interac- 
tions of the force-modification field in the presence of 
non-linear matter fluctuations. It is not sufficient to sim- 
ply have g(a, k) run with wavenumber since the restora- 
tion works in real space and is dependent on the lo- 
cal environment. The PPF approach in |10] for 2-point 
statistics takes the rms density field as an interpola- 
tion parameter between linearized modified gravity and 
non-linear ordinary gravity inside collapsed dark matter 
halos. This approach works quite well for DGP where 
the non-linearity is linked closely with density [37J . The 
original PPF approach requires generalization for other 
non-linear mechanisms which interpolate between the 
linear and deeply non-linear regime in different ways and 
for statistics beyond 2-point functions 

While the non-linear PPF description is still incom- 
plete, it is crucial that a parameterized scheme incorpo- 
rate at least one degree of freedom which allows gravity 
to return to normal locally. Without this feature one 
might infer overly tight limits on modifications to grav- 
ity in the linear regime from non-linear statistics such as 
the cluster abundance or cosmic shear (cf. [38]). This 



4 camb . inf o/ppf 



caveat applies to consistency parameters like 7 as well. 
We turn now to the guidance the non-linear dynamics 
of the worked examples provide on these and related 
issues. 

4. SUPERPOSITION 

Given that the return to ordinary gravity locally ne- 
cessitates some non-linear process, the superposition 
principle for forces among bodies no longer applies. How 
then do we test gravity with cosmological observables in 
the non-linear regime? Does the lack of a superposition 
principle, which implies mode coupling, invalidate even 
large-scale predictions based on linear theory? 

Non-linear Field Dynamics: The most critical and dif- 
ficult piece of a viable description for modified gravity 
is how non-linear effects return ordinary gravity locally. 
The two worked examples f(R) and DGP provide some 
guidance on how this can happen. In both cases, on 
small scales and assuming non-relativistic velocities for 
the matter, Eq. [7] for the lensing potential continues to 
apply but the non-linear equivalent to g comes from a 
modified Poisson equation 

V 2 * = 4irGa 2 p m A rn + ~ V 2 cp , (8) 

where the extra scalar = —df/dR in f(R) theories 
and is associated with the brane position in DGP. Note 
that the Poisson equation is still linear in the density 
fluctuation from the mean A m but gains an extra source 
from a scalar degree of freedom V 2 </>. 

Given that the matter still moves in the metric as 
usual, non-relativistic particles still feel forces according 
to VvP. Nonetheless compared with ordinary gravity 
there is an enhanced gravitational (or "fifth" ) force from 
V(j>. This extra source is set in linear theory to be (f> = 
$ + ^ = g(<I> — "f) and so should obey the equation 

V 2 = g(a)a 2 (8nGp m A m ), (linear, k » aH) . (9) 

To suppress the enhanced forces locally, the non-linear 
generalization should have the form 

V 2 <t> = g{a)a 2 {SnGp m A m - N[4>}) , (10) 

where g(a) is now taken to be the small scale limit of 
the linear relationship between the metric potentials. 
Here N[4>] is some non- linear function of the field and 
its derivatives which goes to 

N[(f>] -> SirGp m A m (11) 

locally. In the f(R) case, N[c/>] — SR(f ((/))) and is a func- 
tion of the local field alone. Here, the general relativistic 
expectation that R = &nGp m gives the minimum of the 
effective potential for <f>. If the minimum can be reached 
inside an overdense region then force law deviations are 
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shielded in the interior by the so-called chameleon mech- 
anism (e.g. [39]). 
In the DGP case 

N[d>] = T l [(VW - (V,V^) 2 ] . (12) 

Here the non-linearity contains quadratic combinations 
of second derivatives of the field just like the Laplace 
operator on the lhs of Eq. (flO)l . Detailed balance sug- 
gests that once V 2 ^ becomes large, the non- linear piece 
will be driven to cancel 87rGp m A m with the field being a 
nearly algebraic function of the overdensity. In this case, 
the suppression is called the Vainshtein mechanism. 

N-body Techniques: The generic features of the non- 
linear dynamics from the worked examples are: the 
Poisson equation for ^ remains linear in its sources of 
the density fluctuation and the Laplacian of the force- 
modification field with the matter falling in the metric as 
usual. Once '5 is obtained other aspects of mesh-based 
iV-body simulations remains the same. 

On the other hand, the scalar source, like the den- 
sity field, obeys a non-linear equation of motion and its 
configuration in the presence of the density field must be 
solved numerically. The non-linear field equation can be 
solved with iterative relaxation techniques [40] . These 
have been successfully applied to the f(R) [44] and DGP 
sources of non-linearity [42 in a particle-mesh N-body. 
Inside a large overdensity, the scalar source responds 
non-linearly and suppresses the additional contribution 
to gravitational forces. Non-linearity implies that the 
gravitational field is no longer the sum of all of the in- 
dependent contributions from individual bodies, instead 
there is a collective saturation of the sources of deviation 
from ordinary gravity. 

This saturation makes scaling relations between lin- 
ear theory predictions and non-linear observables estab- 
lished for ordinary gravity dangerous to apply to mod- 
ified gravity. These include relationships between the 
linear growth rate and the non-linear power spectrum 
and halo mass function [41143142] , 

It is not even guaranteed that the average scalar 
source follows the linear theory prediction on large sup- 
posedly linear scales. Lack of a superposition principle 
breaks the linear assumption of averaging: small scale 
non-linear density fluctuations can affect the behavior 
of the large-scale field. In a universe where even large 
scale density fluctuations are composed from the corre- 
lations between non-linear collapsed dark matter halos, 
this feature can invalidate the predictions of linear the- 
ory. Likewise, the non-linear contributions to the scalar 
field around a collapsed body can alter the field contri- 
bution to the metric substantially. As a consequence it 
can no longer be thought of as moving as a test body 
in an external field. We must solve self-consistently for 



the internal and external field. 

Superposition and Saturation: There are important 
similarities and differences between the ways /(i?) and 
DGP resolve the superposition and saturation issues. 
In both cases, JV-body results show that large scale cos- 
mological density fluctuations do obey linear dynamics 
despite the lack of a superposition principle. 

In the f(R) case, saturation of the scalar source occurs 
as the field itself saturates in value inside an overdense 
region, <fi 0. More specifically, field saturation occurs 
when the gravitational potential of an object exceeds 
the value of the background field </>. 

A thought example where all of the matter resides in 
objects with such deep gravitational potentials exposes 
some counterintuitive properties of the lack of superpo- 
sition. Since the field contribution is suppressed then 
in all regions containing the matter, their contribution 
to the total gravitational potential in Eq. (|8|) reflects 
ordinary gravity. Thus the large scale gravitational po- 
tential will not follow the predictions of linear theory 
with modified gravity. Likewise, the motion of these 
collapsed objects would also not follow the large-scale 
external gradients of the field (set say by the cluster- 
ing of voids). The saturation of the field internal to 
the object would eliminate the influence of any external 
gradient for particles within the object (see Fig. [I]). 

This effect can be phrased as an apparent equivalence 
principle violation between saturated and unsaturated 
objects as pointed out by [H]. That a saturated object 
requires a non-linear field response means that by def- 
inition saturated objects considered as a whole are not 
test particles. Note that the matter that composes these 
objects still are test particles and so these still respond 
to the total metric including the non-superimposable in- 
ternal and external field contributions. As such, there 
is no fundamental equivalence principle violation since 
matter, regardless of its composition, will always fall on 
geodesies of the total, local metric. 

In the cosmological f(R) simulations for cases with 
force law deviations in the linear regime, the dark mat- 
ter halos where most of the matter resides do not possess 
deep enough gravitational potentials to be fully satu- 
rated. In such cases, the linear theory predictions are 
valid and apparent equivalence principle violations be- 
tween dark matter halos of different mass arc difficult 
to identify. 

In the DGP case, the non-linearity involves second 
derivatives of the field. Let us first consider the case 
where the the ViVj</> and V 2 </> terms are comparable, 
and for simplicity proportional [45 , 

A[0] = 4 S ( V W- (13) 



6 



ext 



external 
gradient 



int 



saturated 
field 



+ 



tot 



saturated 
gradient 



Figure 1. Non-linear saturation of field gradients in f(R) (ex- 
aggerated for effect). For a screened object, the internal field <fo nt 
generated by the object does not superimpose with an external 
field 

</>ext- I n particular, the field inside the object loses knowl- 
edge of any exterior gradient and the enhanced gravitational force 
it implies. Screened objects are by definition not test bodies with 
respect to the field but still respond gravitationally to the total 
field </>tot and hence the metric. 



For example, the interior field of an isolated spherically 
symmetric tophat density configuration would have <f> = 
cir 2 + c 2 and hence s = 2/3 [46145] . 

The non-linear field equation then becomes an alge- 
braic relationship between the scalar source and the den- 
sity with the solution J35] 



2gr\ 



y/1 + 32irGg 2 r 2 c sp m A m - 1 



(14) 



At low density V 2 = g(a)a 2 (8nGp m A m ) as desired 
whereas at high density the source is suppressed V 2 oc 

1 /2 

Am . Given the DGP form for g, the threshold for the 
change in behavior scales as A m m 1/s and so the sup- 
pression of the field source is determined by the non- 
linearity of the local density contrast [37] . Contrast this 
with f(R) where the chameleon transition is linked to 
the local gravitational potential. 



In the modified Poisson equation ([8]), this looks like 
a density dependent Newton constant. If this were the 
whole story, DGP would exhibit much stronger viola- 
tions of superposition than f(R). Since much of the 
total mass of the universe is in small collapsed objects 
with A m > 1, these sources would contribute less to the 
large scale modifications to gravity than their average 
would imply as in our thought example above. This is 
problematic for iV-body simulations based on this ap- 
proximation [35] . Even the large scale ( "linear" ) predic- 
tions depend on the resolution: as more and more of the 
mass is resolved into collapsed objects the modification 
to forces on large scales will be reduced. 

Fortunately, this is not the case for the full DGP non- 
linear term. Returning to the isolated source case, the 
field configuration far away from the source (beyond the 
Vainshtein radius) reveals an unscreened source due to 
the ViVj</> terms [45J44] . Hence as long as the linear 
regime scales are larger than the Vainshtein radius of the 
collapsed objects, the lack of a superposition principle 
does not violate linear theory predictions. This fact is 
also related to the absence of an apparent equivalence 
principle violation in DGP [44] . 

Finally, it is interesting to note that for a planar con- 
figuration N[4>] = for DGP and there is no non-linear 
suppression even when A m 3> 1. This has the interest- 
ing consequence that cosmological structure that is pla- 
nar or sheetlike exhibits the largest modification to grav- 
ity. These counterterms that suppress the non-linearity 
for planar configurations are a generic feature of models 
with "galileon" shift symmetry [35] , 

5. DISCUSSION 

Perhaps the primary lesson of the worked examples 
is that gravity is fragile and preserving what we know 
about gravity already poses severe restrictions on a vi- 
able modification that explains cosmic acceleration. On 
the other hand, this fragility provides form to paramet- 
ric approaches and insight on how to design cosmological 
tests that best complement our knowledge of gravity. 

Demanding a metric theory where the matter obeys a 
conservation law requires that deviations take the form 
of a dark energy component under ordinary gravity. The 
most important parameter for modified gravity in this 
language is an effective dark energy anisotropic stress 
that is coupled directly to the lensing potential. The 
resulting phenomenology is therefore distinct from typ- 
ical physical dark energy models. It is important to 
note that the large scale Friedmann and small scale 
Newtonian dynamics implied by the modifications dif- 
fer. They do so for the same reason that dark energy 
cannot be considered smooth relative to the dark mat- 
ter on horizon scales and above: conservation of energy- 
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momentum. The parameterization given in 10J spans 
both regimes and allows linear theory tools for general 
relativity to be used for modified gravity [33 . 

There does not yet exist a complete parameterization 
of the non- linear phenomenology of modified gravity but 
any viable description must contain a mechanism to sup- 
press the modification locally. Density based approaches 
work well for the DGP model but less well for the po- 
tential driven f(R) model. More importantly, the DGP 
and f(R) models provide specific examples of this sup- 
pression. In both cases, the modification is mediated by 
a scalar that obeys a non-linear field equation coupled to 
density fluctuations that can be solved numerically via 
relaxation techniques [40]. Non-linearities in the field 
equation prevent superposition of field solutions given 
the matter. In particular they allow saturation effects 
where the modifications are suppressed within the high 
overdensities where local tests of gravity are performed. 

The drawback of the lack of superposition is that for 
cosmological tests, the field equations must be solved 
jointly with the matter in TV-body simulations. In prin- 
ciple this caveat applies even to large scales in the sup- 
posedly linear regime. The DGP and f(R) models how- 
ever provide examples where the deviations from linear 
theory predictions there are small. 

Likewise, the motion of matter within bound objects, 
while still responding gravitationally only to the met- 
ric, cannot always be thought of as a superposition of 
external and internal influences. These features caution 
against the use of simple scalings that take the linear be- 
havior of modified gravity and seek to predict non- linear 
cosmological observables. 
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